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ABSTRACT 

We use N-body simulations to explore the possibility that the Sagittarius (Sgr) dwarf 
galaxy was originally a late-type, rotating disc galaxy, rather than a non-rotating, 
pressure-supported dwarf spheroidal galaxy, as previously thought. We find that bi- 
furcations in the leading tail of the Sgr stream, similar to those detected by the SDSS 
survey, naturally arise in models where the Sgr disc is misaligned with respect to the 
orbital plane. Moreover, we show that the internal rotation of the progenitor may 
strongly alter the location of the leading tail projected on the sky, and thus affect 
the constraints on the shape of the Milky Way dark matter halo that may be derived 
from modelling the Sgr stream. Our models provide a clear, easily-tested prediction: 
although tidal mass stripping removes a large fraction of the original angular momen- 
tum in the progenitor dwarf galaxy, the remnant core should still rotate with a velocity 
amplitude ~ 20 kms -1 that could be readily detected in future, wide-field kinematic 
surveys of the Sgr dwarf. 
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1 INTRODUCTION 

In spite of extensive theoretical efforts undertaken to re- 
produce the characteristics of the Sagittarius (Sgr) stream 
(e.g. Ibata et al. 2001; Martmez-Delgado et al. 2004; Helmi 
2004; Law et al. 2005, 2010; Fellhauer et al. 2006), there is 
presently no theoretical model that fully explains the wealth 
of available observational data. Two aspects of the stream 
have proved particularly challenging to model. First, the 
leading tail of the Sgr stream is bifurcated, with both arms 
exhibiting similar distances, velocities and metallicity distri- 
butions (Yanny et al. 2009; Niederste-Ostholt et al. 2010), 
which would appear to refute an earlier model explaining 
the bifurcation as two wraps of different ages (Fellhauer et 
al. 2006), or as independent streams from different progeni- 
tors. Second, the position on the sky of the stream suggests 
that the Milky Way (MW) dark matter halo interior to the 
stream has an oblate or perhaps a nearly spherical shape 
(Johnston et al. 2005; Fellhauer et al. 2006), while the helio- 
centric velocities of stream members support a prolate shape 



* Email: jorpega,vasily,nwe@ast. cam. ac.uk 
© 0000 RAS 



(Helmi 2004). In a recent work, Law & Majewski (2010) 
(hereafter LM10) have shown that what appears as mutu- 
ally exclusive results may in fact signal the possibility of the 
halo being triaxial in shape, as Cold Dark Matter (CDM) 
cosmological models predict (e.g. Kazantzidis et al.2010 and 
references therein). However, LM10 find that in order to re- 
produce the location and velocities of the stream the inter- 
mediate axis of the dark matter halo should be aligned with 
the spin vector of the MW disc. This is hard to understand 
as circular orbits about the intermediate axis are unstable 
(e.g. Adams et al. 2007), raising questions as to the forma- 
tion of the Galactic disc. Furthermore, this model does not 
attempt to explain the origin of the stream bifurcation. 

Interestingly, the chemical composition of Sgr clearly 
stands out from the rest of MW dwarf spheroidal 
galaxies (dSphs). First, strong metallicity gradients 
(—2.3,$ [Fe/H] $0.0), unsual in other dwarfs, have been re- 
ported in the remnant core (Giuffrida et al.2010), which also 
extends to the tidal tails (Chou et al. 2007; Monaco et al. 
2007; Keller et al. 2010). Also, recent analysis of the chem- 
ical abundances in the core and tails of Sgr suggest that 
this galaxy underwent an enrichment history more akin to 



Penarrubia et al. 



LMC than to other dSphs (Sbordone et al. 2007; Chou et 
al. 2010). 

Motivated by the above results, here we explore the 
possibility that the Sgr dwarf was originally a late-type, ro- 
tating disc galaxy, rather than a pressure-supported dSph. 
The goal of this paper is to examine how rotation affects 
the properties of the associated tidal stream, and to provide 
model predictions that yield unambiguous tests for this sce- 
nario, rather than reproducing all of the well-documented 
properties of the Sgr stream in detail. 



2 NUMERICAL MODELLING 

The Galaxy model: The MW disc is assumed to follow a 
Miyamoto-Nagai (1975) model with a mass Md = 7.5 x 
1O 1O M0, and radial and vertical scale lengths a = 3.5 kpc 
and b = 0.3 kpc. The MW bulge follows a Hernquist (1990) 
profile with a mass Mi = 1.3 x 1O 1O M0 and a scale radius 
c = 1.2 kpc. The MW dark matter halo is modelled as a 
NFW profile with a virial mass M V ir = 10 12 A/q, virial ra- 
dius r v ir = 258 kpc and concentration c v i r = 12 (Klypin et 
al. 2002). Following LM10, we assume that the halo is tri- 
axial in shape by introducing elliptical coordinates with the 
substitution r — > ra, where m 2 = x 2 /a 2 + y 2 /b 2 + z 2 /c 2 , 
and (a, b, c) are dimensionless quantities. Following the re- 
sults of LM10, we choose the density axis-ratios to be 
(a, b, c) ~ (0.61, 1.34, 1.22). The best-fitting models of LM10 
also suggest that the Sun does not sit on one of the principal 
axis of the halo, but rather slightly off the x (i.e. minor)-axis, 
(xq, 2/q ) = R@(— cos[A], sin[A]), where A ~ 10°. We hold the 
halo parameters fixed through the evolution of our Sgr mod- 
els for simplicity given that tidal streams are barely sensi- 
tive to the past evolution of the host potential (Penarrubia 
et al. 2006). Our analysis also neglects the effects of dynam- 
ical friction on the orbit of Sgr. However, dynamical friction 
is unlikely to introduce a strong orbital decay during the 
time-scale of interest (i.e. the last 2-3 Gyr, see Fig. 5 of 
Penarrubia et al. 2006) 

The orbit of Sgr: The current position of Sgr and its 
line-of-sight velocity in the Galactic standard of rest (GSR) 
frame are (D,l,b) = (25kpc, 5°. 6, -14°. 2) and v los = 171 
km s _1 , respectively (Ibata et al. 1997). To derive the orbit 
of Sgr in the above potential, we adopt the proper motion 
estimates of Dinescu et al. (2005), who find (fii cos b, fib) — 
(-2.35 ±0.20, -2.07 ±0.20) mas/yr. Using the standard so- 
lar motion about the MW centre from Binney & Merrifield 
(1998), this translates into a space motion of (u,v,w) ~ 
(221, —74, 203) km s _1 . We use test particles to integrate its 
orbit back in time for a few orbital periods in order to pro- 
vide initial conditions for our N-body realizations of the Sgr 
dwarf. We find that its orbital peri- and apocentres in the 
above potential are 15 and 67 kpc, respectively, which im- 
plies an orbital period of 1.04 Gyr. Currently, Sgr has just 
(sa 34 Myr ago) gone through its last pericentric passage, 
and will cross the Milky Way disc in approximately 28 Myr. 

N-body realizations of Sgr: We use the algorithm 
buildgal to construct N-body realizations of late-type spi- 
rals composed of a baryonic disc and a dark matter halo (see 
Hernquist 1993 for a detailed description). The disc profile 
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where md is the disc mass, zq is the vertical thickness, and 
Rd the radial scale-length. In our models we assume that 
zo = 0.2Rd- We adopt a non-singular isothermal profile for 
the dark matter halo 
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where rrih is the halo mass, r c is the core radius and a ~ 
1.156. We find that the properties of the Sgr stream model 
are not particularly sensitive to the value of the core radius. 
Here we choose r c = 0.5Rd- 

We use the results of Nierderste-Ostholt et al. (2010) 
to crudely derive a fiducial mass for each of the compo- 
nents of our Sgr N-body models. These authors estimate 
that prior to stellar stripping Sgr had a total luminosity 
of L ~ lO 8 I/0, of which only 30-50% remains currently 
bound in the remnant core. In order to convert N-body 
masses into luminosities we adopt a constant stellar mass- 
to-light ratio of T* = 3.5 (note, however that the strong 
metallicity gradients measured throughout the Sgr core may 
indicate a large range of T*). Under this choice the ini- 
tial Sgr disc mass thus is nid ~ 3.5 x 1O 8 M0. Adopt- 
ing a mass-to- light ratio of T = rrih/L ~ 24, typical for 
dwarf galaxies with L ~ lO 8 !/© (Mateo 1998), we have 
rrih = 2.4 x 1O 9 M0. The total initial mass of our Sgr model 
therefore is M = m d + m h = (T* ± T)L = 2.8 x 1O 9 M . 
Also, since we only simulate the most recent history of the 
Sgr dwarf, our initial conditions must account for the fact 
that the outer halo envelope may have already been lost to 
tides at the time when the stellar stream begins to form (see 
Penarrubia et al. 2008). To do this we impose a truncation 
in the dark matter density profile at r = r cut = 6Rd, which 
roughly corresponds to the tidal radius of a satellite galaxy 
with mass ~ 3 x 1O 9 M0 at a pericentre r pcr i = 15 kpc. 

The only remaining free parameter is the initial scale- 
length of our Sgr N-body models, Rd- We fix its value by 
demanding the bound remnants of the Sgr model to contain 
~ 50% of its initial stellar mass at the final snap-shot of 
the simulation. Indeed, under this definition the value of Rd 
depends on the number of orbital periods for which we follow 
the evolution of the Sgr dwarf. Adopting an integration time 
of 2.5 orbital periods (see §3]), this condition yields Rd — 
0.9 kpc. Under this choice of model parameters, the peak 
rotation velocity of Sgr is 43 km s _1 at R — 2.6Rd — 2.34 
kpc, which are indeed typical values for late-type spirals with 
luminosities L ~ 10 8 Lq (e.g. Swaters et al. 2009). 

The N-body code: We follow the evolution of the Sgr 
N-body model in the Galaxy potential using Superbox, a 
highly efficient particle-mesh gravity code (see Fellhauer et 
al. 2000 for details). Superbox uses three nested grid zones 
centered on the highest-density particle cell of the dwarf. 
Each grid has 128 3 cubic cells: (i) the inner grid has a spac- 
ing of dx — 3-Rd/126 ~ 2.4 x 10 -2 Rd and is meant to re- 
solve the innermost region of the satellite galaxy, (ii) The 
middle grid extends to cover the whole dwarf, with spacing 
lO-Rd/126. (iii) The outermost grid extends out to 50 x r v i r 
and is meant to follow particles that are stripped from the 
dwarf. We choose a constant time-step At = 0.9 Myr, which 
leads to a total energy conservation better than 1% after 
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Figure 1. Leading tail particles in the "Field of Streams" area for different orientations of the Sgr disc with respect to its orbital plane 
(see text). SDSS data are shown in the upper-left panel. Our colour-coding denotes different heliocentric distances (blue, green and red 
for D < 25, 15-40 and 30-60 kpc, respectively). Dashed lines show a projection of the Sgr orbit integrated forward in time. 



the dwarf models are evolved for a Hubble time in isolation. 
This shows that the evolution of the N-body model is free 
from artefacts induced by finite spatial and time resolutions. 



3 RESULTS 

The angle subtended by the internal angular momentum 
vector of Sgr and that of its orbit about the MW, i.e. 
0o = acos[Ji n f Jorb] is a free parameter in our study. The ori- 
entation of the Sgr disc is set at the start of the simulation 
by rotating J; nt an angle 8q about the instantaneous line 
of the orbit. Under this definition, 80 = denotes a model 
where the spin vector of the Sgr disc and the normal vector 
of the orbital plane are perfectly aligned, whilst #0 > and 
60 < respectively indicate models where Sgr rotates in a 
prograde and retrograde motion with respect to its Galactic 
orbit. 

Fig. [l] shows the projection on the "Field of streams" 
area of the sky (Belokurov et al. 2006) of Sgr stellar debris 
for different disc orientations after integrating our N-body 
models for 2.5 orbital periods. For simplicity, we show only 
particles that belong to the leading tail of the Sgr stream, 
given that the trailing tail has not been detected yet in the 
Northern Galactic Hemisphere. This Figure illustrates a few 
interesting points. The first is during the stripping process a 
fraction of the internal angular momentum of Sgr transfers 
to the stream. As a result, the stream tail(s) do not trace 
the progenitor's orbit (dashed lines in Fig. [T]). In practical 



terms, our results suggest that internal rotation in the pro- 
genitor of the Sgr stream may have to be taken into account 
when inferring the shape of the MW halo from fitting the 
position and velocity of stream pieces. A second interest- 
ing point is that bifurcations in the leading tail of the Sgr 
stream naturally appear in this area of the sky if the spin 
vector and the normal vector of the orbital plane are mis- 
aligned. The separation between the bifurcated arms clearly 
becomes more prominent as the value of |#o| increases. In- 
terestingly, by colour-coding the stellar particles according 
to their heliocentric distances, we can appreciate that in all 
models both arms show a very similar gradient throughout 
the sky, in concordance with observational data (Belokurov 
et al. 2006; Niederste-Ostholt et al. 2010). Hence, internal 
rotation in the Sgr dwarf mainly affects the apparent prece- 
sion of the stream plane on the sky. 

The model with #0 = —20° has a striking resemblance 
with the Field of Streams. It also provides a reasonable 
match to most of the existing observational constraints. 
Fig. [3] shows the projected location (upper panel), heliocen- 
tric distance (middle panel) and line-of-sight velocity (lower 
panel) of the leading (open symbols) and trailing (closed 
symbols) tail particles of a Sgr stream model with initial ori- 
entation do = —20°. By colour-coding the particles accord- 
ing to the time at which they are stripped, we can appreciate 
that the bifurcation arises from material lost at consecu- 
tive pericentric passages. In particular, the southern, more 
prominent tail of the leading arm (stream A) corresponds to 
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Figure 2. Projected location (upper panel), heliocentric distance (middle panel) and line-of-sight velocity (lower panel) in the GSR 
frame of the leading/trailing (light/dark colours) tail particles of a Sgr stream model with initial orientation 8q = —20°. Yellow stars and 
open squares denote RR Lyr (Vivas et al. 2005, Watkins et al. 2009) and SDSS detections of Sgr (Belokurov et al. 2006). White particles 
show the remnant core of the Sgr dwarf (see also Fig pi. The dashed line marks the position of the Galactic plane. For ease of reference 
we also show its orbit integrated forward (solid lines) and backward (dotted lines) in time. Note that the bifurcation of the leading arm 
arises from material stripped at different pericentric passages: blue at the penultimate pericentre passage (~ 1 Gyr ago) and red at the 
ante-penultimate one (~ 2 Gyr ago). 



stars that were lost at the third most recent pericentric inter- 
action, i.e. ~ 2 Gyr ago (coloured in red), whilst the fainter 
northern tail (stream B) is more recent and is composed of 
stars that became unbound at the penultimate pericentre, 
i.e. ~ 1 Gyr ago (in blue). This result clearly implies that 
the minimum age of the Sgr stream is two orbital periods, 
which justifies the choice of integration time for our N-body 
models. 

This model correctly reproduces the distances and ve- 
locities measured along streams A and B, as well as the 
recent detections of the trailing tail in the southern hemi- 
sphere. Puzzlingly, although also predicted by our model, 
the presence of the trailing tail in the Field of Streams has 
thus far eluded detection. This may be explained by the 
drop of stars expected beyond the apocentre of the trailing 
tail, i.e. R.A. <; 90°, as well as the sharp increase in distance 
at R.A. ,$120°, which may both conspire to put the sur- 
face brightness and the turn-off apparent magnitude of the 



stream beyond the detection threshold of current photomet- 
ric surveys. 

Fig. [3] suggests that further clues on the nature of the 
Sgr dwarf may be gained by studying its remnant core. In 
the upper panel, we show the projection onto the sky of the 
Sgr dwarf model shown in Fig. [2] Particles are colour-coded 
according to the dwarf surface brightness at their location. 
The core morphology agrees well with that derived from the 
2MASS survey (see e.g. Fig. 4 of Majewski et al. 2003), al- 
though its central surface brightness is somewhat brighter 
than estimated by these authors (/10 ~ 24.8mag/" 2 ). Mis- 
matches in surface brightness, however, may be a conse- 
quence of adopting a constant stellar mass-to-light ratio for 
our N-body particles. The lower panel of Fig. [3] shows that, 
although tidal mass stripping removes a large fraction of the 
original angular momentum in the progenitor Sgr disc, the 
remnant core is predicted to rotate with a velocity amplitude 
~ 20 km s _1 , which translates into a net velocity difference 
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indirectly affect any constraint on the shape of the MW halo 
derived from pressure-supported dwarf models. 

Fortunately, before we embark upon more complex 
modelling of the Sgr stream, there is a clear-cut way to test 
whether the Sgr dwarf was indeed a rotating galaxy. Al- 
though tidal stripping efficiently removes angular momen- 
tum from the progenitor dwarf, we find that the remnant 
core should still rotate with a velocity amplitude close to 
~ 20 kms -1 given the current estimates of the fraction of 
light residing in the tidal stream. Validating this prediction 
is feasible with existing instruments and will shed light on 
the true origin of the Sgr dwarf, as well as on the shape of 
the MW dark matter halo. 

JP and MNO acknowledge financial support from the 
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Kingdom, whilst VB thanks for the Royal Society for fi- 
nancial support. This work has benefited greatly from dis- 
cussions with D. Lynden-Bell. We thank the referee, D. Law, 
for his thorough read of this manuscript and his useful com- 
ments. 



Figure 3. The remnant core of the Sgr dwarf model shown in 
Fig. [2] Particles are colour-coded according to the dwarf surface 
brightness at their location (upper panel) and mean line-of-sight 
velocity (lower panel). Symbols denote the location of the fields 
spcctroscopically surveyed by Ibata et al. (1997). Note that our 
models predict that the remnant core still rotates with a velocity 
amplitude ~ 20 km s — 1 . 



of ~ 40 km s _1 between stars at declinations above and be- 
low 8 ~ —30°. Although not shown here, the central line 
of sight velocity dispersion is ~ 12 km s _1 , in reasonable 
agreement with the estimates of Bellazzini et al. (2008), who 
find a uniformly flat profile at a ~ 10 knr s _1 over the central 
0' < r sC 9' range. 

The only available kinematic survey throughout the Sgr 
dwarf dates back to fbata et al. (1997). These authors mea- 
sured a line-of-sight velocity between 160-170 km s" 1 in 
most of their pencil-beam fields, hence with little or no evi- 
dence of net rotation. Although in these locations our model 
accurately reproduces their velocity measurements, Fig. [3] 
suggests that Ibata et al. fields were too sparsely sampled 
to pick up the rotation signal predicted by our models. 



4 CONCLUSIONS 

Here, we explored the possibility that the Sgr was originally 
a late-type, rotating disc galaxy. We show that bifurcations 
in the leading tail of the Sgr stream, similar to those de- 
tected by the SDSS survey, naturally arise in models where 
the disc is misaligned with respect to the orbital plane. This 
occurs because material is primarily stripped at pericentric 
passages, and successive passages occur at different orien- 
tations of the Sgr disc. Together with observed metallic- 
ity patterns that are atypical for dwarf spheroidal galaxies, 
this suggests that the Sgr dwarf may have originally been a 
galaxy akin to late-type spirals with a peak rotation of c 45 
kms - . We also find that internal rotation alters the posi- 
tion of the stream with respect the Sgr orbit, which nright 
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